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We reconsider the geometry of pure and mixed states in a finite quantum system. The ranges 
of eigenvalues of the density matrices delimit a regular symplex (Hypertetrahedron Tjv) in any 
dimension A'^; the polytope isometry group is the symmetric group Sjv+i, and splits Tjv in chambers, 
the orbits of the states under the projective group PU{N + 1). The type of states correlates with 
the vertices, edges, faces, etc. of the polytope, with the vertices making up a base of orthogonal 
pure states. The entropy function as a measure of the purity of these states is also easily calculable; 
we draw and consider some isentropic surfaces. The Casimir invariants acquire then also a more 
transparent interpretation. 

PACS numbers: 03.65.Ca, 02.10.De, 02.40.Et 



I. INTRODUCTION 

Pure states tt in Quantum Mechanics are described by 
rays in a Hilbert space Ti., or equivalently as unidimen- 
sional subspaces. As such, they can be identified with 
orthogonal projectors of unit trace (so for fc-dimensional 
subspaces they will have trace fc); denoting by End V the 
hncar operators of a vector space V, we have 

{space of pure states} w {tt e End H\ n ~ n'' , 

Tt[tt] = 1, TT^ = vr} (1) 

To include mixed states p, it is enough (von Neumann) 
to relax idempotency ( tt^ = tt, so Spectrum tt — {1,0}) 
to mere positivity, p > <-!■ Spec p > 0, which impUes 
< Pi < 1 for Pi G Spec p. Hence 

{space of pure/mixed states} ~ {p G End Ti\ p = p^ , 

Tr[p] = l, p>0} (2) 

So the pure states are extremal, as p^ < p, and one 
can write the mixed states as incoherent addition of pure 
ones. 



^PiTT,, ^Pj = 1, > 



(3) 



For finite quantum systems { H ^ V, dim V = N < 
+00 ) the geometry of these states has been studied lately 
from different points of view, as it is pertinent to quan- 
tum information and other modern appUcations. For 
example, in [l| the nature of the space of pure states, 
namely the projective spaces CP^~\, is stressed, as they 
are rank-one symmetric spaces; in [2| the mixed states 
are related to a coherence vector ~n, lying in a subset 
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of euclidean space, determined by the values of Casimir 
invariants for the group SU{N). Q studies the stratifi- 
cation of general states under the projective group 



PU{N) = U{N)/U{1) = SU{N)/Zn 



(4) 



in relation to a (in general not regular) foliation by Kahler 
manifolds, similar to the orbits of the (co-)adjoint repre- 
sentation of SU{N), that is, the Kirillov-Souriau method. 
The geometry of the N = 3 case had been studied some- 
time ago also by Michel [3] in relation to SU{3) as a flavor 
group. 

In this report we would like to complete the above de- 
scriptions in several ways. In particular, we elaborate 
an idea in j^, showing that the set of eigenvalues pi of 
the density operator p, p^ G Spec p, coincides precisely 
with the points of a solid simplex (the hypertetrahedron 
Tn), the simplest regular polytope. The natural isome- 
try group, the full symmetric group Sn+i, corresponds, 
of course, to a (finite) Coxeter group, that is, the Weyl 
group of the Lie algebra of the group SU{N+1): the 
unit trace in our case is traded for the traceless character 
of (anti-)hermitian matrices describing A^. The orbit 
space Tn/Sn+i is stih an (irregular) hypertetrahedron, 
of size 1/{N+1)1 of the previous regular one, behaves like 
the Weyl chambers, and describes precisely the orbits of 
our states under the projective group, each point being 
just an orbit. The type of orbits, as classified by the lit- 
tle groups, is in correspondence with the combinatorial 
elements of the simplices, namely vertices, edges, faces 
and so on, and also with the partitions of the number N. 

For example for iV = 2 the set of pure states is the 2- 
sphere = CP^ , the set of all states is the three dimen- 
sional Bloch ball B'^ of unit radius, with boundary being 
the pure states, the simplex is just the closed segment I: 
[-1, 1], the Weyl group Z2 is reflection in the middle, and 
the chamber is the half-segment I/Z2 =[0, 1] with the 
most mixed states (called mix Max or mM henceforth), 
and with 1 indicating the (sphere of) pure states. In 
this simplest (and nonrepresentative) case there is only 
one stratum, i.e. one type of orbit (namely 2-spheres), 
besides of course the fix point mM or 0; the (single) 
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Casimir I2 just labels the radius of these spheres Sr , 
< r < 1. The mM state corresponds to partition [2], 
the rest to [1^], as stratum orbits « [/(2)/f/(l)2. 

The Casimir invariants h admit a double interpreta- 
tion, as an homogeneous system of generators of the cen- 
ter of the enveloping algebra of Lie group, or in our case, 
as the symmetric functions (Newton) over the roots of 
the spectral equation for the density operator. We elab- 
orate here the considerations of [J] , based on the pioneer 
work of Biedenharn Q . 

In this picture it is also clear how to compute the en- 
tropy rj of these mixed states, where 7y(p) = - Trplog(p). 
This varies between 77 = log N for the mM state to log 
1=0 for pure states. We include some graphs showing 
the entropy for some boundary lines on the "chambers" 
of the orbit space, as well as some isentropic surfaces in 
the general case. Each case N includes, in a precise sense, 
all the previous ones, n < N. 

We stress also the action of the projective isometry 
group PU{N + 1) C SO{{N + l)^ - 1) as the group 
acting effectively in the space of states, in the sense 
of the characterization of geometries (F. Klein): in the 
transformation p — > UpW the effective group acting is 
PU = U/U{1): the kernel of the U action is t7(l), and 
even SU acts with kernel Z„: the explanative diagram is 





-> SU{N) - 


PSU{N) 


I 


i 


II 


U{1) 


^ UiN) - 


PU{N) 


I 


I 




U{1) 


== U{1) 





The organization of the paper is as follows. In Sect. 2 we 
establish notation and show several properties of our ob- 
jects needed later. In Sect. 3 we recall the situation for 
qubits (A^=2), qutris (iV=3), and N—A, incorporating 
the new observations of above and suggesting general- 
izations to higher dimensions. Then Sect. 4 deals with 
the general theory for arbitrary N: we exhibit explicitely 
the polytope, the Weyl quotient ("chamber"), the types 
of states, the Casimir invariants and discuss the entropy 
function, and also some cases of isentropic surfaces. 

Other considerations (including mention of omissions) 
are in our final Sect. 5. 



II. GENERAL DESCRIPTION AND 
PROPERTIES OF MIXED STATES 

In this Section we sumarize important results of to 
[3|)- Let S be set of TV-dimensional density matrices p, 
and let P C 5 the subset of pure states, tt, with tt^ = n. 
The conditions 

p = p\ Tr[p] = l (6) 

are respectively linear and afhne, so the whole set of 
solutions of the equation is like a real vector space of 



dim N'^-1. The positivity condition p > 0, though, is 
a convex condition, and selects a nonlinear submanifold 
S of (which is closed as the spectral restriction 

is really < pi < 1), of the same dimension, with V 
as the extremals of the convex set (not the boundary, in 
general: we know of course that V = CP^^^, so dim 
V = 2N - 2 <s: - 1 for large N). In the set of her- 
mitian traceless matrices {h} one introduces the definite 
scalar product 

</ii,/i2>:-(l/2)Tr[/ii/i2] (7) 

and therefore there exists an orthonormal base Ai,A2, 
Xn2_i with properties 

\ = Xl Tr[A,]=0, Tr[A2]=2. (8) 
So the general density matrix can be written as 

where the factor y/N{N — l)/2 guarantees that pure 
states, TT^ = TT have norm \n\ — 1; see [3] for details; 
n is called the coherence vector. 

We recall also two related properties of the group 
SU{N); first, the square of the adjoint representation 
contains the adjoint BOTH in the symmetric and in the 
antisymmetric part (for N^2 the first is missing), and 
second, consequently, this induces two algebras in the 
space of matrices (corresponding to the coefficients fijk 
and dijk of Gell-Mann; of course, the /'s correspond to 
the Lie algebra structure); for details see Q. 

Now the group PU{N) is acting in the set of states 
iS as p UpU^^. If the eigenvalues of p are differ- 
ent, the little group is U{1)^~^, and then the generic 
orbit is PU{N)/U{1)^^\ of dimension iV(iV-l). This 
is a Kahler manifold, as are all orbits of the adjoint of 
any simple Lie group (we notice that N{N-1) is even); in 
particular this space is called a (complex) flag manifold, 
with structure (as homogeneous space), 

¥lc{N) = PU{N)/U{1)^-' = U{N)/U{1)^. (10) 

The other extreme contains the pure states, with little 
group U{N— 1); these most critical orbits have dimension 
2(iV-l), as we said above. There is always naturally a 
single most-mixed state, O (or mM), with n = 0: this 
is the fixed point under the PU action, and it is unique. 
Orbits in state space with conjugate little groups form 
what is called a stratum (see e.g. [4|), and in fact for a 
compact group (as is our case) acting "nicely" in some 
space X the number of strata is finite; we shall see that 
it coincides with the number of types of states and also 
naturally with the partitions of N. 

We define the entropy rj = r]{p) of a state p as the 
expectation value of the operator H— -log [p) > (von 
Neumann; notice p < 1, so log(p) is negative); recalling 



3 



Pi mM P, mM P 

I 1 1 I 1 

FIG. 1: "(color online)" Eigenvalue segment for N — 2 

that in the density formahsm < A >= Tr(pA), we have, 
with /ij G Spec p, 

r^{p) = -Trplog(p) = \og{IipT^^^). (H) 

This varies, as said, between log N and log 1=0: the mM 
state is the most disordered, and the pure states are the 
most ordered: we shall see also the N case reproduces 
the entropy function for all previous cases n < N, and 
therefore there must be several manifolds of mixed states 
isentropic, with the same value for the entropy: we shall 
provide some examples. 

As for the Casimir invariants, we define them from the 
spectral equation for the density operator as the coeffi- 
cients of the powers, i. e. : 

- (Trp)p^-i + /2p^-2 - /3P^-^ + ...{-Ifdetp = 

(12) 

where Ii = Tr (p) = 1, (the first Casimir). They deter- 
mine the spectrum up to a permutation, e.g. 
-^2 = J2i<j i^ii^j^ a-nd /„ = Il/ii = Hi...p.N 
These operators can be expressed also as traces of pow- 
ers of representative matrices, as for example 

h = (l/2)[(Trp)2 - Tr(p2)]; so it's zero iff p is pure. 
h - (l/6)[(Trp)3 + 2Tr(p3) - 3Tr(p)Tr(p2)]. 

There are also several inequalities assuring all Casimirs 
are nonnegative, etc. 

III. THE SIMPLEST CASES 

For N=l we have just a single point, a pure state, 
as CP° = point, with entropy log(l)= 0. For N=2 the 
density matrix can be written as 

p=(l/2)(l + a.x); (13) 

we have yi € ball (radius 1) with the 2-dim bound- 
ary sphere — CP^ of pure states. The set of eigen- 
values {(1 + z)/2,(l — z)/2} or {x,y) with x + y = 1 
makes up a segment I in fig [H where Pi is e.g. the 
state [0, -1], mM is (1/2, 1/2), and P2 is [1, 0]. The 
crucial but trivial point now is that under the symme- 
try Z2: {x,y) — > (y,a;) the two states are equivalent, so 
the segment I becomes just the half-segment [0, 1]; that 
is, the mM state (0), the mixed states (0 < x < 1) 
and the (representative of) pure state(s) F, x— 1. In 
this case PU{2) = SU(2)/Z2 = 50(3), so the orbits 
are just spheres, and x in the half-segment means just 
the radius of them. There are only two strata, the 
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FIG. 2: "(color online)" Entropy for A'^=2 case 



fixed point (little group SO{2>)) and the rest, little group 
S0(2)=U(1)). The two strata correspond to the two 
partitions namely U{2)/U{2) = {Point} ^mM is [2] 
and C/(2)/[/(l)2 = PU{2)/U{1) = CP^ is the partition 

The entropy is a smooth function from log(2)=0.693 
(for the mM state) to log(l)=0 for the pure one (fig 
[2]): we shall see that the entropy function of the N case 
always contains that of the n < N previous cases. 

There is only a Casimir, the quadratic one, with 

h = [(1 + z)/2][(l - z)/2] = xy^ x{l ~ x), (14) 

which lies between 1/4 and zero. 

One can also use an "angle" picture [2^, namely with 
z = cos9, < 9 < n, X = cos^e/2,y = 8111^9/2. The 
(quadratic) Casimir is now (l/4)sin^0 . The entropy in 
terms of angle variables is 

77(61) = -cos2(6l/2)log[cos2(6l/2)]-sin2(6'/2)log[sin2(6l/2)]. 

(15) 

For N=3 we have p = (1/3)(1 + V3 n • A* ), where the 
A's are e.g. the Gell-Mann matrices. Choosing A3 and Ag 
diagonal, we get, with n = (0, 0, a, 0, 0, 0, 0, b) and 

I f 1 + \/3a + b \ 

pMag = l~Via + b (16) 

V l-2b J 

s^s pdiag = (l + \/3(aA3-|-6A8))/3. Hence positivity implies 
— 1 < b < 1/2, and one gets a regular triangle "upside 
down" inscribed in the circle of radius one (fig [3|) . 

Notice the three vertices A, 2, 3; 3 edges ^1,2,3, and 
the three interior lines (heights) /ii,2,3 with intersection 
points Qi.2,3. The figure is identical to the so-called 
Fano plane, related (among other things) to the octo- 
nion multiplication rule [6] and to the projective plane 
F2P^,with automorphisms the simple group PS'L3(2)of 
168 elements [7]. This is another test of the relation of 
SU{3) with octonions! 

Under the action of the symmetry group, now S3 with 
6 elements, isomorphic to the dihedral group D3 (the 
isometry of the regular n-gon is D^), the fundamental 
domain or chamber is a little rectangular triangle (there 



FIG. 3; "(color online)" Positivity domain for 7V=3 
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are six of them), 1/6 in area of the big one; we select one 
of them as AOQP as in fig. [3l Here O is the mix Max 
state, P the pure, Q intermediate, i. e., r?(0)=log(3), 
7y(Q)=log(2), 77(P)= log(l)=0. 

In terms of the eigenvalues of p, O is (1/3, 1/3, 1/3), 
P is (1, 0, 0), Q is (1/2, 1/2, 0); recall that the order is 
immaterial. The line OQ is type h {xxy, x >y); the line 
QP is type / {xyO, x > y > 0), and this line reproduces 
the N=2 case, with Q seen as the mM state. 

The line OP is type h, with {xyy, x > y). Note all the 
statements are order-invariant! In fact in the large trian- 
gle the line POQ is a single line! Notice also the bound- 
ary lines in the big triangle unite pure states, so they 
have in the spectrum. Each line represents the three 
embeddings of SU{2) in SU{3), called by Gell-Mann u, v 
and w spin. It becomes a single line in the little triangle, 
line PQ. 

Note the set of generic states (little group U{1)'^) is 6- 
dimensional, open and dense; therefore, the representa- 
tives are bi-dimensional (2-|-6=8=dim SU{3)). The crit- 
ical states (two equal eigenvalues, little group U{2)) are 
the two lines emanating from O, but O is the fixed point, 
with isotropy just PU (3) itself. The three partitions of [3] 
are: mM = U{3)/U{3) is [3], U{3)/U{1) x U{2) = CP^ 
and U{3)/U{1)^, which is the flag manifold. Notice also 
the little triangle is rectangular, NOT regular. 

The Casimir invariants now are two, quadratic and cu- 
bic. They are 

^2 = AiiM2 + Ai2Ai3 + AisMi; h ^ ^^l^^2^J■3 (17) 

The quadratic Casimir is zero for the pure state, the cubic 
zero for pure state also and on the boundary line QP. We 
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There is an inequality between I2 and /a, because the 
cubic equation has to have real roots; see [J|. 



FIG. 4: " (color online)" The entropy surface over the triangle 
OPQ 




o 

FIG. 5: " (color online)" Constant entropy lines over the tri- 
angle OPQ 

One can also parameterize the states with two angles 
with e.g. 

X = sin^ 61/ 2 cos2(/)/ 2, ?/ = sin^ 61/2 sin^ 0/2, z^cos'^e/2 

(18) 

The entropy of the (xyz) state, with x + y + z = 1, 
is \og{x^^y^'^ z~^). Fig [4] shows the entropy surface. It 
varies from log(3) at O to log(2) at Q to log(l)=0 at P. 
Fig [5] shows the isentropic lines over the triangle. The 
point R on line OP with ?y(i?) = log(2) is xyy =.768, 
.116, .116. 

We have come now to the N = 4 case. The general 
state with Ai, . . . A15 is 

p-i(l + %/6n.A). (19) 

We take A15 as 

Ai5 = 4Fdiag{l,l,l,-3}, (20) 
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which gives the diagonal density matrix 



/I 









1 - V6a + V2b - 









1 - 2V2b 








1 - 3c/ 



(21) 




FIG. 6: "(color online)" The little tetrahedron for N=4: 



in particular c < 1/3. The figure is now a regular tetra- 



hedron inscribed in S^. Again, the action of the 24 el- 
ements of the 84 group generates the little rectangular 
tetrahedron as quotient (figd]). 

In the eigenvalue notation (xyzu), x + y-\-z-\-u~ 1, 
the coordinates of the vertices are: 

O is (1/4, 1/4, 1/4, 1/4); mM statc=C/(4)/C/(4), par- 
tition [4]. 

P is (1, 0, 0, 0); pure state=C/(4)/[C/(l) x C/(3)], par- 
tition [3,1]. 

Qa is (1/2, 1/2, 0, 0); nongeneric mixed 
state=J7(4)/C/(2)2, partition [2,2]. 

Qf is (1/3, 1/3, 1/3, 0); nongeneric mixed 
state=C/(4)/[J7(l) x [/(3)], partition [3,1]. 

There are six lines joining vertices, as follows 





OP 


OQa 


OQf 


QaP 


QaQf 


QfP 


Length 
Spectrum 
type 


1 

X > yyy 
[3,1] 


i/Vs 

XX > yy 
[2,2] 


1/3 
XXX > y 
[3,1] 


VV^ 

xyOO 
[2,1,1] 


V2/3 

XX > yO 
[2,1,1] 


2V2/3 

yy > xO 

[2,1,1] 



There arc four faces, w e describe just two here: 





OPQa 


PQfQa 


Spectrum 
Type 


xyzz 
[2,1,1] 


xyzO 
[1,1,1,1] 



Finally, there are five types of orbits: 



Spectrum 


1/4, 1/4, 1/4, 1/4 




xxyy 


xyzz 


xyzu 




Representative 


mM=0 


y=iooo 


Qa 


QaP 


Interior 




dim 





6 


8 


10 


12 




character 


Fixed point 


Pure states « CP^ 






Generic (Flag manifold) 




G/H 


U(4)/U(4) 


U(4)/[U(l)xU(3)] 


U(4)/U(2)2 


U(4)/[U(2)xU(l)2] 


U(4)/U(l)4 





The discussion of the Casimir invariants and the calcu- 
lation of the entropy follows similar lines to the = 3 
case, so we omit the details here. The face QfQaP re- 
produces the previous N = i case. This N — A case 
is important as being the first case with correlation, for 
example entanglement. 



IV. GENERALIZATION 

The density matrix for an N level system is given by 

where the A's are A^^ — 1 Hermitian traceless N x N ma- 
trices with square 2; A^ — 1 of them can be diagonalized 
simultaneously; let us call them A3, Ag, A15, AAr2_]^. 
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Then the density matrix is given by 



1 

N 



I l + V^ (a + VV3 + c/V6 + ... + zy=3=:) 














.) 







l-(iV-l)z / 
M 



So again, —1 < z < 1/{N — 1); the ranges of eigenvalues 
always give us a regular simplex or solid hypertetrahe- 
dron T/v, because there are + 1 orthogonal pure states 
in CP^ , symmetrically distributed in the S^^^ sphere of 
unit radius. The hyperfaces of the polytope delimit the 
range. The center corresponds to the mM state, which 
remains invariant under U{N + 1) or PU{N + 1). 

The (full) symmetry group is Sn+i, of course (arbi- 
trary permutation of vertices etc.); this divides T/v in A'^! 
rectangular, irregular little hypertetrahedra. In any of 
them we have O, the mM state; P, the vertex; and A^ — 1 
Q's, from Qp (center of the hyperface) to Qa (center of 
the edge uniting, in Ta?, P with another vertex V'). 

Notice in this picture some elements are interior in the 
original, regular polytope, others are at the boundary 
(extremal); for example, the cell which does not contain 
the O — mM is in the boundary, but the others are 
interior; similar for other elements: faces, edges, etc. 

The calculation of the Casimir invariants starting from 
eq. (23) is mechanical and trivial. Some general results 
are 

In{— detp) ~ for a boundary state (24) 
because of a zero eigenvalue. Also 

In-1 = for a boundary edge state (25) 
because of two zero eigenvalues, etc. Besides 



Ij{mM state) = 




because all the eigenvalues are equal to 1/N. 

Let us discuss the entropy in this general situation; the 
formula 

77(p)=log(n/i-'^-), ^/i. = l (27) 

can be applied without any difhculty; again, we just in- 
clude some results 

77(0) = logiV > r){p) > log(P) = log 1 = (28) 

for O ^ p ^ P. Besides for A= edge, F= face, we have 

7^{Qa) = log2, 7]iQF) = log3, . . . , riiQceii) = log(A^-l) 

(29) 



where Qa is in the boundary edge, in the boundary 
face, Qceii in the boundary cell, the {N — l)-polytope. 
There are plenty of isentropic surfaces, which we refrain 
to make explicit, as we already established them in the 
iV = 3 case. See also [l]. 

V. FINAL REMARKS 

Our purpose in this work has been to describe the 
eigenvalue set and the orbit space of density matrices 
in a finite quantum system. The picture is 

{set of states} solid hypertetrahedron Tjv (30) 
{orbit space} = T^/Sn+i ^ rectangular 

nonregular small /i-tetrahedron tj^ (31) 

where t]\r = {set of density states} /PU{N). 

The actual geometry of mixed states, out of the N — 
1,2 cases, is rather involved [8] and we have not pre- 
tended to improve on it. There is also the interesting 
problem of finding maps between density states, posi- 
tive/strictly positive, untouched also here; see e.g. the 
recent work '9'|. 

The geometry of the complex flag manifolds 
U {N)/U{1)'^ is very rich; we remark here only one result 

U{3)/U(lf := Flc{3) = SU{3)/U{lf - CP^ CP^ 

(32) 

where the twisted product is similar to the used in [lO| 
to express the homology of Lie groups. 

The N = 3 case reminds one of the 3-dim Jordan alge- 
bras: the action of SU{3) on the 3x3 hermitian trace- 
less matrices is similar to the action of the exceptional 
group P4 on the exceptional Jordan octonionic algebra; 
indeed, the Moufang plane OP^ is F4/Spin{9) (Borel, 
1950), whereas our CP^ is SU{3)/U{2). In fact, the 
two cases are related, as the following beautiful chain 
of groups show 

SU{3) Spin{8) Eq 

a/ dt a/ n a/ (33) 
SO{i) G2 F4, 



7 



where a =Aut leads to the fixes sets (i.e. Eq F4, C. (Spain) for grant FPA-2006-20315. 
Spin{8) G2 and SU{3) SO{3) ), and t (triali- 
tiy) and d (duaUty) means to relate Spin(8) with and 
SU{3) with G2. That is, F4 « Spin{8) plus the 3 S-dim 
representations and G2 ~ SU{3) plus 3 and 3. 
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